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ABSTRACT 


H. optimal control theory, based on singular value loop shaping, ts used to 
synthesize a controller for the statically unstable longitudinal dynamics of X-29 
aircraft. Two design cases are studied: 2-input 2-output; and 3-input 3-output cases. 
H. theory provides a direct, effective procedure for synthesizing control laws satisfying 
specified performance objectives and robustness specifications. The 2 I/O case has 
better performance, a faster response and is more robust, than the 3 [T/O case. 


Discussion and comparison of results are given. 
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I. INTRODUCTION 

The main task for control system designers is to synthesize a control law which 
maintains system response and error signals to within prespecified tolerances while at 
the same time ensure relatively high robustness. Some measures of system robustness 
commonly used are stability margins, sensitivity and disturbance attenuation. Thus, a 
robust system denotes a system with satisfied system performance, has large stability 
margins, good disturbance rejection and low sensitivity to parameter variations. 

For the single-input single-output (SISO) systems, some techniques, 1.e., root- 
locus method, Bode plot and Nyquist diagram, had been well developed to design and 
estimate system performance and robustness. Some of these classical methods can be 
generalized to multi-input multi-output (MIMO) systems, however, not all of them do. 
For instance, from a Nyquist diagram, one can not get a satisfactory notion of 
multivariable stability margins. 

In the past decade, a lot of research on properties of multivariable feedback 
system has been accomplished. Singular value Bode plots of return difference and 
loop gain matrices have emerged as useful indicators of multivariable robustness [Refs. 
1, 2, 3]. Among the techniques that have been recently developed, H., frequency- 
weighted LQG, and LQG loop transfer recovery optimal synthesis neones have made 
singular value loop shaping a routine matter [Ref. 5]. For a multivariable feedback 
control system, singular value loop shaping involves the manipulation of system loop 
gain over a specified frequency bandwidth to tmprove system performance and 


robustness. 


H,, theory provides a direct, reliable procedure for synthesizing a controller. 
While the frequency-weighted LQG optimal synthesis theory (or "H, theory”) and LQG 
loop transfer recovery theory lead to somewhat less direct but highly effective iterative 
procedures for shaping singular value Bode plots to satisfy singular value loop shaping 
specification [Ref. 5]. 

The aim of this thesis is to present a controller design procedure, bases on the 
H.. theory, for a reduced order, linearized longitudinal dynamics model of the X-29 
aircraft analog backup mode. Some review about the properties of singular value, 
sensitivity reduction and return difference matrices will be given in Chapter II. In 
Chapter I, H,, theory and weighting matrix are introduced. Application of H, theory 


to X-29 longitudinal dynamics model is made in Chapter IV. The design result is 


discussed in Chapter V. 


fl. FEEDBACK PROPERTIES 


For a control system, some properties such as impulse response matrix, transfer 
function matrix and other characterizations of the response to commands can be 
changed by prefiltering the command signal, rather than using feedback loop. But 
some properties, usually termed feedback properties, like stability, sensitivity and 
disturbance attenuation can be altered only through the feedback [Ref. 2]. Here, we 


focus on these feedback properties which also define the robustness of system. 


A. RETURN DIFFERENCE MATRIX 

Consider the standard linear time-invariant multivariable feedback control system 
in Fig. 2.1. It consists of the interconnected plant (G), controller (F) and sensor (H) 
forced by commands (Fr), measurement noise (n), and plant disturbance (d). The 
prefilter (P) is an optional element used to achieve deliberate command shaping. 


From this configuration, we get 


plant 

y(s) = G(s)u(s) + d(s) 221) 
controller 

u(s) = P(s)r(s) - F(s)z(s) (2.2) 
sensor 

z(s) = H(s)(y(s) + n(s)) (23) 


where y(s), d(s), m(s), u(s), z(s). r(s) are vectors. and G(s), F(s), H(s) are matrices. 
For the feedback loop, define the transfer matrix 


M(s) = F(s)H(s) 
3 





Figure 2.1 Standard Feedback Control System 


and matrices 
L,(s) = G(s)M(s) 
L,(s) = M(s)G(s) 
are called return ratio matrices [Ref. 2] at the y(s) mode and the u(s) node, 
respectively. 
And 
I + L,(s) 
I + L,(s) 
are return difference matrices at y(s) node and u(s) node. Simuarly, the y(s) node and 


u(s) node inverse return difference matrices are defined, respectively, as 


Fone. 

I+ L,’ 
Those terms introduced above are fundamental and important in the control system 
synthesis. Because for a SISO system, the phase and gain margins can be determined 
by the behavior of return difference as a function of frequency. And it tums out, as 
we will see later, the return difference matrix also provides a means of measuring 


robustness in a multivariable system. 


B. SINGULAR VALUES AND MATRIX NORM 
In the design of a feedback loop, one needs to determine the "size" of matrix 
to estimate its properties. For a vector X(x,, X,, X3,.-.,X,), the size can be defined as 
the Euclidean norm 
WX ova? 4x2 4x7 4 000 + x) 

For a matrix, a more general concept is needed. One way to describe the size of 
matrix A is based on the largest singular value of A, as 

6(A) = max( Il AX Il / Il X Il ) 
where Il « Il denotes Euclidean norms. The singular values 6, of a rank n matrix Ae 
C™ are defined as the non-negative square roots of the eigenvalues of A°A, where A’ 
is the transpose of the complex conjugate of A. It is convenient to order them as 
follow 

Gre. 2.0..... 2c Gp 


If rank r < n, then there are nm - r zeros singular values. that is 


One uses singular values rather then eigenvalues to describe matrix size since 
eigenvalues can be a very poor indicators of the "size" of matrix [Refs. 5, 14]. The 
singular value is also a good measure of the near-singularity of a matrix, i.e., how 
near the determinant of a matrix comes to being zero. The condition number, the ratio 
of the largest to the smallest singular value, provides information about sensitivity to 
perturbation, which also is an important indicator of robustness. Listed as follows are 
some useful properties of singular value : 

(1) G(A) = max (Il AX Il / Il X Il 

(2) g(A) = min (ll AX Il / Il X Il) (the least singular value of A) 

(3) o(A) Ss | A(A) | S GCA) , where A, denotes the i-th eigenvalue of A 
(4) If A” exists, o(A) = 1 / G(A") 

(5) If A” exists, o(A) = 1 / o(A") 

(6) o(AB) Ss o(A) o(B) 

(7) 2%°.., 07 = Trace (A‘A) (Ref. 5] 

In SISO system, it is customary to define sensitivity in terms of percentage 
variation in system output resulting from a given percentage change in the plant with 
the input held fixed. A similar way is used to develop the sensitivity function for a 


MIMO system. 


C. SENSITIVITY FUNCTION AND DISTURBANCE ATTENUATION 
Consider a finite dimensional linear time invariant (FDLTI) system in Fig. 2.2, 

where r, d, and n are command, disturbance, and measure noise respectively. If it is 

asymptotically stable, then we have input-output loop error and sensitivity relation as 


follows : 





Figure 2.2 Finite Dimentional Linear time Invariant System 


input-output 
y = GF(I + GF)'(r - n) + (I + GF)'d (2.4) 
loop error 
e=r-y 
= (1 + GF)'(r - d) + GF( + GF)'n (2.5) 


system sensitivity 
AH, = (1 + G'F)'AHo. (2.6) 
where AH, and AH), represent the changes in the close-loop system and nominal 
open loop system caused by changes AG in the plant G. ie, G’ = G + AG [Ref. 1]. 
Here, we define 


S(s) = (I + L(s))" = e/r (27) 


R(s) = F(s)(I + L(s))' = uff (2.8) 

T(s) = Lis) + L(s))’ = I - S(s) = yr (2.9) 
where L(s) = G(s)F(s) is the loop transfer function, S(s) 1s known as the sensitivity 
function matrix, T(s) is complementary sensitivity function. It can be seen that loop 
transfer function L(s) determines the other three matrices. And, as will be seen, the 
singular value bode plots of these matrices play important role in the design of 
multivariable robust control system. 

From Eq. (2.5), it is seen that the loop error " e " resulting from disturbance 
inputs can be made small by making the sensitivity matrix, or inverse of the return 
difference matrix (I + GF)', small. From Eq. (2.6) the system sensitivity to plant 
variations may be also reduced by decreasing the sensitivity function. One can 
decrease the size of the sensitivity function by making 

o(SG@)) = o[(1 + LGo))"] 
small, or (from singular value properties (5)) 

oI + Lg@)) (2.10) 
large. 
Note also 

o(I + LGw)) = a(Lgo)) for G(LGw)) » 1 
This tells us that large loop gains or " tight " loops yield good performance. But one 
can not make loop gains arbitrarily high over arbitrarily large frequency. Since we can 
see, from Eq. (2.5), large o(L(jm)) over a large frequency range not only make errors 


due to r and d small, but make error due to n large. Large gains can also lead to 


instability. There are certain performance tradeoffs between command / disturbance 


error reduction and sensor noise error reduction and stability to be satisfied. 


D. STABILITY AND STABILITY MARGIN 
In SISO system, gain margin and phase margin are the stability margins that are 
commonly used. Both are conveniently expressed in terms of the magnitude of the 
retum difference, I + LGw). In MIMO system, the retum difference retains its 
importance, but the concept of gain margin and phase margin become problematic. 
The stability of the system in Fig. 2.2 is directly related to the matrices S(s), 
R(s) and T(s). To find the stability margin for the system is to determine a lower 


bound on the " size " of smallest perturbation to the plant that will destabilize the 
system. And the singular value Bode plot of S(s) and T(s) can be used as the measure 


of stability margins. 


i. H. CONTROL DESIGN 
H_ control theory provides a direct, reliable procedure for synthesizing a feedback 
controller designed to meet singular value loop shaping requirements. Its capability 
includes addressing the full range of stability margins, sensitivity, and robustness 
optimization and response-shaping problems that can be formulated within the singular 
value Bode plot framework. The standard configuration of an H, problem is shown 
in Fig. 3.1. The design objective is to find a stabilizing controller F(s) for the 


augmented plant P(s) while keeping the robustness specifications satisfied. 


P(s) yi 


F(s) 


Figure 3,1 Standard H. Small Gain Problem 


10 


The purpose of this chapter is to present a brief overview of H_ optimal theory 


and a control system design procedure based on this theory. 


A. CONTROL SYSTEM FORMULATION 

As shown by Eq. (2.5), (2.6), and (2.10), the loop error due to disturbance and 
system Sensitivity to plant vanations can be reduced by suppressing the sensitivity 
function S(s). (There is also some stability requirements desired in the robust system 
synthesizing.) These goals can be achieved by properly selecting the weighting 
functions in the design procedure. 

Shown in Fig. 3.2, 1s the compensation configuration of system in Fig. 2.2, where 
P(s) is the augmented plant with weighting functions W,, W, and W, which penalizing 


the error signal, control signal and output signal respectively. 


Augmented Plant P(s) 


y2 


Figure 3.2 Compensated System with Augmented Plant P(s) 


11 


From this Figure, we get 


Wie W, -W.G ql 
yi 1 Viv 0 W, u, 
y2 | ic 0 W,G 


af) tft Se | 


Where the augmented plant with weighting functions is denoted as 


W, -W,G 
0 W, 
P(s) = 
0 W,G 
hi? @ake iq 


A state space realization of P(s) is given by 


A B, B, 
P(s) = C, D,, D., (3.1) 
C, D,, D,. 


(3.2) 





| 73 


If the weighting transfer function is improper (i.e., has more zeros than poles), as can 
be the case for W,, then no state space realization exists. The produce W,G is in 
this case still proper and can be represented by a state space realization. Some 
calculations are need to compute the state space realization of the augmented plant. 
From Fig. 3.1 one can define the transfer function T,,,, as 
¥, = Ty at, . 
This transfer function can be expressed in terms of the weighting matrix and previous 


defined loop matrices, i.e., 


Vii Ws 
Yiv = W.R | u, 
Yi- W,T 


And the close loop transfer function from input u, to output y, is denoted as 
WS 
w= W.R 
W,T 
Where S, R, and T are loop matrices defined in Eqs. (2.7), (2.8),(2.9). 
The compensated system is said to be internally stable if the augmented A matrix 
of the compensated system is stable, i.e., when external input signal u, equals to zero, 


the states of both P(s) and F(s) will go asymptotically to zero for any initial 


conditions. And the controller is said to be stabilizing. 


Is 


B. H, OPTIMAL CONTROL PROBLEM 
The H,, optimal control problem is defined as the following small gain 
problem: 
For the given transfer function matrix P(s), find a stabilizing controller F(s) 
such that the close loop transfer function matrix T,,,, is internally stable and its 
infinity-norm is less than or equal to 1, Le., 
! 


ll T _ <a (3.9) 


ylul 
Where the H.. - norm of matrix A is defined as 
| A il, = sup G(AGa)). (sup : the least upper bound) 


" 


This problem is called a " small gain problem " because the closed loop gain is small 
(1e., less than or equal to one). As mentioned before, the largest singular values of 
S(s) determine the plant disturbance attenuation, since S(s) is the same as the close 
loop transfer function from disturbance d to plant output y (Fig. 2.2). One may 
specify the system disturbance attenuation performance with the frequency-dependent 
weighting function W,, that is 
o(S(jw)) < | W,"(Ga) | 

Now, consider the system in Fig. 3.4, where A,(s) and A,, represent the additive 

plant perturbations and the multiplicative plant perturbations respectively. Defining the 


size " of A,G@) and A,(j@) as G(A,(j@)) and GOA Gm), the stability theorems are 


obtained as follows: 


PERTURBED PLANT 





Figure 3.3. Additive / Multiplicative uncertainty 


Robustness Theorem I: 
If the system in Fig. 3.3 is stable as both A, and Ay are zero. Let A, = 
0, the size of the smallest stable A,{s) for which the system become unstable 
1s 
J 
G(AGwo)) = ———— 
o(TGo)) 
Where T(jw) is the complementary sensitivity function of the system. It is seen that 
decreasing the value of o(T(jm)) will increase the size of A(jq@), that is, increase the 


stability margin. 


Robustness Theorem 2 : 


15 


Suppose the system in Fig. 3.3 is Stable when A, and A,, are both zero. Let 
Ay = 0. Then the size of the smallest stable A,(s) which destabilizes the 
system 1s 
_ ] 
o(A,Go)) = ———_—_—- 
o(RY@)) [Ref. 5] 
We can see that the smaller o6(R(j@)) is, the greater will be the stability margin. 
With properly selected weighting function matrices W,(jw) and W,(o), the 
stability margins of the control system can be specified by the following singular value 
inequalities, 
o(R(jm)) < | W,'Ga) | 
G(T(jo)) < | W,'(ja) | 
Together with the performance specification W,Gq@), the design specifications of the 


robust control system may be written as 


o(S(jo)) < | W,'Go) | (3.3) 
G(R(ja)) < | W,'Go) | (3.4) 
6(T(jo)) < | W,'Go) |. (3.5) 


Thus, for a robust control system synthesis, the requirements may be presented as 
reducing the plant disturbance effect as much as possible while meeting the control 
inputs and the robustness constraints. This goal can be achieved by the singular value 
loop shaping iterative procedure that will be introduced in the following sections. 

As a consequence of singular value properties and the definition of H, norm, 
these singular value inequalities can be combined into the single infinity norm 


specification in terms of the close loop transfer function T,,,, of the form 


eae, lo sal 


ylul 
where 

WS 
is W.R 
W,T 


Note, the matrix T,,,, is partitioned into submatrices representing performance 


requirements and stability constraints. 


C. WEIGHTING MATRIX SELECTION 
Since the plant disturbance can be attenuated by reducing the value of O(s(ja)) 
continuously, a constant scaler weighting gain y is added to W, to facilitate the 
singular value looping shaping process. That is 
a(S(jw)) < | y' * W,"(ja) | 
So, aS we increase the value of y, the sensitivity function 1s suppressed down 
continuously. On the other hand, from Eq. (2.9) (T(s) + S(s) = I, we know that T(s) 
goes to I as y increases (since S(s) goes to zero). This means the complementary 
Sensitivity is forced against its upper constraint W,'(s) as the value of y increases. 
The closed loop bandwidth, therefore, is widened. This improves the system 
performance, 1.e., make system response faster. At the same time, the singular values 
of closed loop transfer function T,,,, 1s pushed up to upper limit. The singular value 
loop shaping process stops as the condition in eq. (3.9) is reached. 
It is interesting to note that, 
S(s) = (I + L(s))' = L(sy"' 


as O(L(j@)) » 1, ( ie., the system has large loop gains) 
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1 / o(SGo@)) = o(LGo)) 
Similarly, as o(LGj@)) « 1, 

T(s) = Ld + L(s))' = L(s) 
that is, 

o(TGw)) = o(LGo)) 
These can be seen in Fig. 3.4, ie., above the zero dB line, g(L(jw)) » 1, 1 / o(SGo)) 
= 6(L(jm)), while below the zero dB line, 6(LG@)) « 1, o6(TGo@)) = o( LG@)). Thus 
the robust system design consists of the use of high loop gains at low frequencies 
where the performance specifications are defined and the use of low loop gains at high 
frequencies where the robustness constraints lie [Refs. 5, 13]. The system performance 
will be limited within the range between 6(G(s)F(s)) curve and 9(G(s)F(s)) curve. 

The weightings W, and W, are selected as transfer function matrices of which 

diagonal elements are frequency-dependent constants. And the size of these weighting 
matrices is the same as the number of plant output states. This small gain problem 
is solved by using a software program, ie., " hinf.m " provided by MATLAB 
ROBUST-CONTROL TOOLBOX package. The number of states of the stabilizing 
controller F(s) produced by " hinf.m " is the same as that of P(s). In addition it is 
required by " hinf " that the D,, submatrix of P(s) (see Eq. (3.1)) is full column rank. 
An easy way to ensure that this is the case is to choose the weighting matrix W, with 
a invertible " D-matrix ", e.g., W,(s) = eI, where € is any non-zero number [Ref. 5]. 
The parameters considered in selecting these weightings are gain. order, attenuation 


and comer frequency, which are used to achieve small o(S(s)), i.e., large loop gains, 
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Figure 3.4 Singular Value Spec’s on S(s) and T(s) 


over the broadest possible bandwidth subject to the robustness constraints. And it 
should be noted that the controller F(s) depends on the weighting matrices selected. 
Once the weightings are selected, the control system design is in a computer-aided 


design environment and requires only " one parameter " y-iteration. The flow chart 


of the H_ iterative procedure is shown in Fig. 3.5. 
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Figure 3.5 HL Singular Value Loop Shaping Iterative Process 
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IV. H. CONTROLLER SYNTHESIS FOR X-29 FIGHTER 

The purpose of this Chapter is to present a modern aircraft controller design 
procedure based on the H, optimal control theory. The model used is the longitudinal 
dynamic model of X-29 fighter. The X-29 is a single seat forward swept wing (FSW) 
demonstrator aircraft build by Grumman Co., which is designed to be a new generation 
of tactical fighter with the advantages of light weight, low cost and high efficiency. 
The computer-aided design program, MATLAB-ROBUST CONTROL TOOL BOX, is 
used in controller synthesis. The brief plant description and design objective are given 
in the following sections. The design procedure and results are also presented. The 


script files specifically written or modified for this problem are listed in appendix A. 


A. PLANT DESCRIPTION 

The aerodynamic advantages of forward swept wing design include improved 
maneuverability, with virtually spin-proof characteristics, better low-speed handling and 
reduced stalling speed. The FSW X-29 (see Fig. 4.1) was selected primarily because 
its multiple, independently controlled surface make it an ideal candidate for multiloop 
synthesis of advanced control mode. 

The aircraft is designed to have a -35% stability margin at subsonic speed. A 
wing strake extends aft from the trailing-edge at each wing root, each strake with a 
trailing-edge flap which has its own integrated servo actuator to augment foreplanes 
for pitch control. All-moving canard surfaces. one on each side of the center fuselage, 


outboard of engine inlet ducts and operated by servo actuators for primary pitch 
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Canard 


Figure 4.1 FSW X-29 Fighter 


control. In addition the aircraft 1s equipped with flaperons (Ref. 10]. 

The X-29 longitudinal dynamic model considered in this study is an analog 
reversion mode with the aircraft trimmed at 0.5 mark, 30,000 feet. The design point 
was one of the critical nodes in the analysis of the aircraft. The 14-state reduced- 
order model includes a short period approximation of the aircraft longitudinal dynamics, 
vertical velocity w and pitch rate q, and 3 fourth order actuator dynamics for the 
longitudinal control surfaces respectively, i.e., the canards, flaperons and strakes. 

For the first design, a 2-input, 2-output configuration is used. Two separated 
commands r, and r, are input to the three control surface actuators with r, controlling 


the canards and strakes and r, controlling the flaps. The outputs of the system are two 


Pape 


aircraft states, w and q. The physical configuration of this 2-input, 2-output open loop 
actuator/aircraft dynamics model is shown in Fig. 4.2. 

The second design example considered is a 3-input, 3-output controller synthesis. 
Three separated commands r,, r, and r, are input to the three control surfaces canards, 
flaps and strakes respectively. The output to be controlled are w, q and the canard 
control input 6,. The physical configuration for this model is shown in Fig. 4.3. 
Comparison is made between this two proceeding designs. 

The state variables are listed in Table 4.1. The state space realization of the 
open loop plant model G 1s listed in Appendix B. The poles of this system are listed 
in Table 4.2. It is seen that there is a positive pole, ie., 1.9550, on the real axis. 
That means the X-29 has an unstable short period mode. 

The 14-state model was scaled to improve the numerical conditioning of its state 
space representation. The w state was transformed to angle of attack o using the initial 
forward velocity uy, 1e., @ = w /u,. The units of the actuator third derivative states 
were transformed from rad/sec’ to le+04 rad/sec’. This scaling was effective in 
reducing the condition number of the system state matrix A, from an order of 
magnitude of 10'° to 10* [Ref. 13]. 

The singular value plot of the uncompensated plant G(s) is presented in Fig. 4.4. 
The solid curve and dashed curve represent the maximum and minimum singular values 
of G(s) respectively. As can be seen, at low frequency, the open loop plant has poor 
disturbance attenuation, high sensitivity to plant variations and modeling errors, and a 
small control bandwidth. These properties will be improved by singular value loop 


shaping procedure presented in the following sections. 
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Figure 4.2 X-29 2-Input, 2-Output Open Loop Configuration 


: W. L. Rogers [Ref. 13] 
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X-29 3 I/O Open Loop Configuration 
Source 


Figure 4.3 


Table 4.1 Discription of The Uncompensated X-29 state variables 


State Description 
a ang le-of —-attack 
gq pitch rate 
6. canard control input 
é, flap econtrEet. ino ie 
é. strake control input 
é canard control rate 
b, f lacweonerolmmace 
é. strake control rate 
é. canard control accel. 
oy flap control accel. 
é. strake control accel. 
5. Canara controle qenk: 
5, flap ‘controler. 
é. strake control jerk 


Source : W. L. Rogers [Ref. 13] 
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Units 


rad 
rad/sec 
rad 
rad 
rad 
rad/sec 
rad/sec 


rad/sec 


rad/sec* 


rad/sec* 


rad/sec* 


le+04 rad/sec> 


le+04 rad/sec° 


le+04 rad/sec° 


Table 4.2 Uncompensated X-29 open loop poles 


SV - dB 


Figure 4.4 Singular Value Plot of The Uncompensated Plant G(s) 


10-5 


poleg = 


—-2.2746e+02 
-2.2746e+02 
-1.4455e+02 
-1.4491e+02 


-1.0031e+02 
-2.7155e+00 
-~-5.2518e+01 
=5 925 18e701 
-~-5.2506e+01 
-5.2506e+01 
-5.0067e+01 
-2.0172e+01 
-2.0115e+01 


1.9550e+00 


2.3201le+021 
2.3201le+02i 


4.8255e+011 
4.8255e+011 
4.8410e+011 
4.8410e+011 
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B. DESIGN SPECIFICATIONS AND WEIGHTING FUNCTIONS 

Since the design objective is to attenuate the disturbance effect as much as 
possible within loop gain constraint and robustness specifications. The design 
specifications, in this case, are : 

Performance Spec.: 

Minimize the sensitivity function as much as possible. 

Robustness Spec.: 

1). Attenuate the close-loop singular values of the complementary 
sensitivity matrix by 20 db at frequencies beyond @ = 100 rad/sec 
(15.9 Hz) (This will ensure that the system will have sufficient stability 
margin to tolerate variations the loop transfer function of magnitude as 
large as a factor of ten at frequency w = 100 rad/sec). 

2). a second-order roll-off beyond wm = 100 rad/sec (This ensures the 
controller is proper and matches the uncontrolled drop-off). 

With these considerations, we have covered the most important issue in feedback 
control system design. Once the design specifications are decided, the selection of 
weighting matrices is the next step. For the 2-input 2-output system, the weighting 
matrices are chose as follows : 


O1(1+s/ 01 ) 


OWS)! = ¥ 2) ae (4.1) 
(1+s/100 ) (2 x 2) 
W.~(s) = - 0.025 * I. (4.2) 
(4 x 4) 
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1000 
EG) = (4.3) 
oe (2 x 2) 


*? ft 9 


to penalize the error " e ", control input " u " and output " y " respectively. 

The singular value plot of W,'(s), W,'(s) weighting functions is given in Fig. 
4.5. The robustness specification W,'(s) has 20 dB drop-off and 40 dB roll-off at 
w = 100 rad/sec which satisfied the design requirements. As mentioned, W,'(s) and 


W,,'(s) should be so chose to ensure the 0 dB crossover frequency of W, is sufficiently 


lower than the 0 dB crossover frequency of W,7’(s). W,(s) is selected to penalize the 


X-29 Design Specifications 
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Figure 4.5 Singular Value Plot of Weighting Functions 
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control input u, and the weighting factor 0.025 is chosen to confine the actuators’ 
response within their design limits. The selected W,(s) also ensures that the submatrix 
D,, of P(s) has full column rank as required in the two Riccati equation solution to the 
H.. problem. The H, theory requires that D,,(s) matrix is full column rank in order 
for an optimal H,, controller to eliminate infinite energy inputs. If we decrease the 
weighting factor, more energy will input to the system, in this case, the actuators’ 


response may exceed the limits [Ref. 13]. 


C. DESIGN PROCEDURE 

Before starting the singular value loop shaping procedure, it is pec to find 
the balanced state space realization for the unaugmented plant for a better numerical 
condition. This can be seen from the different condition numbers of state matrix A, 


before and after balanced realization, which are 5.526e+04 and 2.5567e+02 respectively. 


1. 2-INPUT, 2-OUTPUT DESIGN CASE 

With the selected weighting matrices, the augmented plant P(s) of the 
balanced plant is of order eighteenth, W,(s) and W, each adding two states to the 
system G(s). Since W,(s) is not a proper transfer function Eq. (4.3) (e., it has more 
zeros than poles), so it has no state space realization. But W,(s)G(s) is proper and can 
be realized in state space form. 

With the computer-aided design program, y is chosen as | for the first try 
of the iteration process. As continuously increasing the value of y, the singular value 
of cost function T,,,, 1s pushed up to its upper limit gradually until no solution exists 
for any larger y. After several iterations, yY equals to 2.6, and the maximum singular 


values of T,,,, reached the all pass limit (1.e., 0 dB line). Fig. 4.6 is the singular value 
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COST FUNCTION Tylul (Gamma = 2.6) 
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Figure 4.6 Cost Function T,,,, for y = 2.6 


plot of cost function T,,,,, the solid line represents the maximum singular values and 
dashed line represents the minimum singular values. Fig. 4.7 and 4.8 are the singular 
value plots of complementary sensitivity and sensitivity, the dashed line and dotted line 
again represent the maximum and minimum singular values respectively. The 
complementary function is pushed flat against the specified limit W,’, as the value 
of gamma is increased, the closed loop bandwidth is widened (to 30 rad/sec) and the 
response speed is increased. Similarly, the sensitivity function is pushed down against 
its limit y'W,' as y value increased. The singular value of I + G(s)F(s) loop gains 
are increased to about 50 dB at low frequency for this design. which markedly 
improves the X-29 performance properties. This large loop gains indicate good 


disturbance attenuation and low sensitivity to plant variations over a control bandwidth 
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COMP. SENSITIVITY FUNCTION AND 1/W¥3 
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Figure 4.7 Complementary Sensitivity Function T(s) and W,' for y = 2.6 
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Figure 4.8 Sensitivity Function S(s) and W," for y = 2.6 


a 


of about 3 rad/sec and satisfy the stability requirements. Note, in Fig. 4.8, W,"(s) = 
o(S(j@)) at low frequency. And there is about 40 dB drop-off from @ = 10 rad/sec to 
@ = 0.1 rad/sec which is the range of interesting. The state space realization of the 
18-state stabilizing controller F(s) is listed in Appendix C. Two uncontrollable states 
are removed using PRO-MATLAB function " minreal.m ". A fast transient poles is 


eliminated while retaining the most important characteristics of the closed loop system. 


Thus, the final controller has 15 states. The poles of this reduced-order controller are 


listed in Table 4.3. 


Table 4.3. Poles of The Reduced-Order Controller 


poleacpr 


-2.1954e+02 

-6.5170e+01 + 
-6.5170e+01 
-1.4487e+02 
-1.3577e+02 


9.5444e+011 
9.5444e+011 


-5.3024e+01 + 4.9271le+01li 
=o. 50/2GC+O1 —24.29207 eto 11 
—-4.1572e+01 + 3.1587e+01i 
-4.1572e+01 - 3.1587e+01i 
-1.3180e+01 + 1.9602e+01i 
-1.3180e+01 - 1.9602e+01i 
-3.0118e+01 + 1.8361le+01i 


-3.0118e+01 
-1.0000e-02 
-1.0000e-02 


so 


1.8361e+011 





Figure 4.9 X-29 Feedback Configuration 


Fig. 4.9 shows the closed loop system of compensated X-29, where Fs) is 
the 15th order controller and G(s) 1s the 14th order X-29 plant. Since F(s) is placed 
in series with plant G(s), the command vector r is composed of the reference 
commands r, and r, to controlled outputs a and gq. So , the closed loop compensated 
system has 2 inputs, 2 outputs and 29 states. The poles of this reduced-order closed 
loop system are listed in Table 4.4. The singular value plot of output, inverse-retum 
difference matrix I + (G(s)F(s))" (also called the output, multiplicative retum difference 
matrix) for the compensated system is shown in Fig. 4.10. The minimum singular 
value of this matrix, ie., 0,1 + (G(j@)F(j@))']. provides a measure of the aircraft's 
gain and phase margins with respect to multiplicative modeling errors using the 


universal gain and phase margin curve (Fig.4.11) (Ref. 4]. 
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Table 4.4 Poles of The 2 I/O Closed Loop System 


—-2.2739e+02 + 2.3194e+02i1 
—-2.2739e+02 - 2.3194e+02i 
-2.1723e+02 
-6.4876e+01 + 9.6310e+01i 
-6.4876e+01 - 9.6310e+01i 
-1.3558e+02 
-1.4595e+02 
-1.4469e+02 + 9.2844e-O02j 


-1.4469e+02 - 9.2844e-02i1 
-1.1118e+02 

—-4.9782e+01 + 4.9587e+01i 
-4.9782e+01 - 4.9587e+01li 
-5.2188e+01 + 4.8038e+01i 
-5.2188e+01 - 4.8038e+01li 
-5.2464e+01 + 4.8363e+01i 
-5.2464e+01 - 4.8363e+01i 
-4.1968e+01 + 3.2086e+01i 
-4.1968e+01 - 3.2086e+01i 
-4.1045e+01 

-2.7932e+01 + 1.3447e+0O1li 
-2.7932e+01 - 1.3447e+0O1i 
-1.1994e+01 + 1.1959e+01i 


-1.1994e+01 
-1.9603e+00 
-2.6697e+00 
-3.8935e+00 
-1.0123e+01 
-1.6232e+01 
-1.9727e+01 


1.1959e+011 


It is seen in Fig. 4.10 that o,,,[1 + (G(s)F(s))’] drops to about -3 dB at 


frequencies between (0.3 rad/sec and 20 rad/sec. Entering the vertical axis of the 
universal curve at singular value of 0.71 (1.e., -3 dB). The gain and phase margins of 
compensated X-29 near 0 dB crossover frequency are -11 dB to +5 dB and +42 deg 
respectively. This is more stable than the -8 dB to +4 dB and +35 deg gain and phase 


margins for the typical design of a fighter aircraft. 
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X-29 SV PLOT (I + inv(GF)) 
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Figure 4.10 Singular Value Plot of The Output Inverse-Retum Difference 
Matrix for 2 I/O System 
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Figure 4.11 Universal Gain and Phase Margin Curve 
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Shown in Fig. 4.12 are the closed loop Bode plots of a and q respectively. 
The solid curves are the response to r, and dashed curves are to r,. It is seen that, at 
low frequency, a is “all pass" (i.e., has 0 dB gain) to input 1 but has very small 
response (-50 dB) to input 2. Similarly, as seen that q is “all pass" to r, but has small 
response to r,. This indicates the decoupling of a from q to r, and r,. Note, the 
System performance at high frequencies satisfy the robustness requirements for both 
cases. 

The time response of a and q to 1 deg/sec (0.01745 rad/sec) step input are 
presented in Fig. 4.13 and 4.14. It is seen from Fig. 4.13 that the response of o 1s fast 
while with rise time of 0.3 seconds and 0.1 second to r, and r, respectively. We also 
can see that the decoupling of m and q to the two commands from both plots. If we 
reduce the weighting function W,(s), more energy will be input to system, and that 
energy increases the response speed and the decoupling will be more pronounced. 

Referring to Fig. 4.15 and 4.16, the time response of control deflections 6,, 
6,, and 6,. It is seen that the peak magnitudes are 0.013 rad, 0.006 rad, and 0.007 rad 
for 6,, 5, and 6, respectively, which are well within the X-29’s control surface 
deflection limits listed below : 

(1). canards (leading edge) : 30 deg up / 60 deg down, 

(2). flaps (trailing edge) : 10 deg up / 25 deg down, 

(3). strakes (trailing edge) : 30 deg up and down. 
The peak control rates (see Fig. 4.17 and 4.18) are 0.0013 rad/sec, 0.0004 rad/sec, and 
0.0008 rad/sec respectively, which are much smaller than the X-29 actuator minimum 
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Figure 4.12 Closed Loop Bode Plot of « and q for 2 I/O System 
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Figure 4.13 Step Response of @ for 2 I/O System 
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Figure 4.14 Step Response of q for 2 I/O System 
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Figure 4.15 Step Response of Canard Control Input 6, and Flap Control 
Input 6, 
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Figure 4.16 Step Response of Strake Control Input 4, 
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Figure 4.17 Step Response of Canard Control Rate 6, 
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Figure 4.18 Step Response of 5, and 3, 
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(1). canards : 1.75 rad/sec 
(2). flaps : 0.87 rad/sec 
(3). strakes : 0.52 rad/sec. 
The response is well within the actuator limits but one also has a very small step- 
input. 
2. 3-INPUT, 3-OUTPUT DESIGN CASE 
The second design example we use is the same X-29 dynamic model as in 
the previous example, but with different input and output vectors, 1.e., r,, r,, and r, to 
control canards, flaps and strakes separately and outputs are a, q, and 6, (see Fig. 
4.3). Thus, it is a 3-input, 3-output system with 14 states. 
The same design specifications are used for this case, and the weighting 
functions are chose as follows : 


O1¢(1+s/.01 ) 





(NSIS) = Vt ee (4.4) 
(1+s/ 100 ) (3 x 3) 
Wes) = -01025 *_L (4.5) 
(6 x 6) 
1000 | 
Wre(s) = * |] (4.6) 
s? (3 x 3) 


The augmented plant P(s) is an 20th order system with W,(s) and W,(s) each adding 
three states to the plant G(s). Following the same iterative procedure as of the 


previous example, the cost function reaches its all pass limit as y = 1.3 (see Fig. 4.19), 
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Figure 4.19 Cost Function T,,,, for y = 1.3 


ylul 


then the iterative process stops. The solid, dashed, and dotted curves in this plot 
represent the maximum, medium, and minimum singular values of T,,,, respectively. 


The controller F(s) produced has 20 states, that is the same as number of states of 


ee t 


P(s). Three states are removed using minreal.m function and one state is 
eliminated by " obalreal._m " program. So, the final controller has 16 states, the 
controller / plant series has 30 states in total (uncompensated plant G(s) has 14 states). 
The poles of this reduced-order closed loop system are listed in Appendix E. 
Comparing Fig. 4.20 with Fig. 4.7 and 4.8 shows that this 3 I/O system are not as 
robust as the 2 I/O model. The 3 I/O system has smaller disturbance attenuation, 


larger sensitivity to plant variations and modeling errors, smaller control bandwidth and 


closed loop bandwidth. Fig. 4.21 is the singular value plot of output, inverse-retum 
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Figure 4.20 Complementary Sensitivity T(s) and Sensitivity S(s) for y = 1.3 
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Figure 4.21 Output, Inverse Retum Difference Matrix for 3 I/O System 


difference matrix. It is seen that the o,,,{ I + ((G(s)F(s))' ] (ie., the dotted curve) 
drops to approximately -6 dB between 0.2 rad/sec and 13 rad/sec. Enter the vertical 
axis of the universal curve (Fig. 4.11) at singular value of 0.5 (i.e., -6 dB), the gain 
and phase margins near 0 dB crossover frequency are -6 dB to +3 db and + 28 deg 
respectively. It indicates that the stability margins of the 3 I/O system, with respect 
to multiplicative modeling errors, are less than the values desired of a fighter aircraft, 
i.e., -8 dB to +4 dB, +35 deg gain and phase margin respectively. Obviously, it is 
also not as stable as the 2 I/O system discussed before. 

From Fig. 4.22, the Bode plots of a and q to the 3 inputs, we also can see 
that a is decoupled from q to r, and r,. And it is interesting to see that both o and 


q have very little response to r,. This also can be seen from Fig. 4.23, 1e., r, has 
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Figure 4.22 Closed Loop Bode Plots of a and q for 3 I/O System 
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Response of Alpha to 0.01745 rad/sec Step Input 
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Figure 4.23 Step Response of o and q for 3 YO System 
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almost no effect on both a and q. This indicates that there is not much difference 
to system response between 2 inputs design or 3 inputs design. 

Shown in Fig. 4.24 and 4.25, the peak values of control deflections 6,, 6, 
and 6, (1.e., 0.005 rad, 0.0017 rad, and 0.015 rad respectively) are all well within the 
design limits listed before, and are also smaller than those of the 2 I/O system since 
the loop gains are smaller. Similar results are obtained from the plots of control rates 
(Fig. 4.26 and 4.27), ie., all responses stay within their design limits. 

Thus, we may conclude that the performance of this 3 I/O control system satisfies 
the design specifications but is not better than those of 2 I/O system designed 


previously. 
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x10-5 Response of DC to 0.01745 rad/sec Step Input 


Canard Input - rad 
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x10-3 Response of DF to 0.01745 rad/sec Step Input 
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Figure 4.24 Step Response of 6, and 6, for 3 I/O System 
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Response of DS to 0.01745 rad/sec Step Input 
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Figure 4.25 Step Response of 6, for 3 I/O System 


x10-4 Response of DCDOT to 0.01745 rad/sec Step Input 
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Figure 4.26 Step Response of 5. for 3 I/O System 
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Figure 4.27 Step Response of 5, and 5, for 3 I/O System 
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V. CONCLUSIONS AND RECOMMENDATIONS 

Although the mathematic work behind the H,, theory is very complicated, the 
implementation of this method is easy using a well developed computer-aided design 
software. From the design procedure presented above, it is seen that H.. method 
provides a direct, systematic and effective way to synthesize controllers for the 
statically unstable longitudinal dynamics of X-29 aircraft. The two-command control 
design example demonstrates a fast response, large stability margins, good disturbance 
attenuation and low sensitivity performance. While the 3-input control configuration 
shows a somewhat less performance due to the smaller loop gains. 

It is seen in both design cases that the two control outputs, angle-of-attack o 
and pitch rate q, are decoupled to the command inputs. This is an unportant feature 
which is needed to effect the advanced control modes in which aircraft flight path and 
pitch attitude are independently controlled. Since, as shown in the 3-input case, the 
separated commanded strake has almost no influence on the control outputs, we may 
conclude that a reasonable and effective controller can be based on the use of two- 
command control in the X-29 control system . 

In both designs, we found the peak values of control surface deflections for the 
actuators and control rates are all much smaller than their design limits. The selected 
W, weighting matrix ensures this compatibility but at a severe penalty in performance. 

Further reduced order modeling of the controller is possible but was not 


accomplished in this thesis. 


a 


APPENDIX A 


Matlab Script Files 


diary nx29.dat 
format short e 


disp(’ e) 

aisp(’ This script file is designed to solve the Hinf optimal control ‘) 
aisp(’ problem for the X-29. The 14th order FDLTI model, in state ’) 
disp(’ space form, is that of the X-29 aircraft and actuator dynamics. ’) 
disp(’ Two states are those of the aircraft dynamics, i.e., alpha ’) 
disp(’ and g. The remaining 12 states are the dynamics of the ’) 

disp(’ three, fourth order actuators, i.e., the canard, flaperon, and ’) 
disp(’ strake actuators. The order of the unbalanced states is as ’) 
disp(’ follows: ’) 

adisp(’ ia) 

disp(’ alpha, q, dc, af, ds, dcdot, didot, dsdetj dedbidoee,) 

disp(’ dfdbldot, dsdbideot, dctrpldot, dftepldot, astrpiace ) 

disp(’ ae) 

disp(’ Given the open loop transfer function G(s)=Cinv(Is-A)B+D, a ’) 
disp(’ stabilizing controller F(s) will be found such that the Hinf norm’) 
disp(’ of Tylul is less than or equal to one. '’) 

disp (’ ‘) 

disp(’ $3 

$pause 

ele 

disp(’ X-29 aircraft and actuator state space representation’) 

disp(' ‘) 

diss) 

ag={~-.4181d+00 -9960a+00 -.2269d-01 -.1213da+00 <-.1948dad-01 ~.9493da-03 


-44270-04 ~-.6712d-04 .14514-05)4-.21626-04 —.35400d=05 5) 0.050). Omron, 
~-5474d+01 -.3424d+00 ~2o000rO] —. 13I0CAO+Oly  —-. 10550701 SOc G mer 


-.1164d-01 -.6397d-02 -.2509d-03 -.5362d-03 -.2912d-03 0.0 0.0 0.0; 
0.0 0.0 0:30. 0.0°960.0 1.0 0.0990.0 §0.050.0 3020) 0.0 OU wece. 
0.0 0.0 0.0 Q20 O70" 0.0 8120 0.0 "0-0 0.0 070 ~ 0-0) GO, Ue. 
0. Ong050. 2040 0.0.50. 0% 30.02 -0. 0 sa 0 m0 Ome 0. One OO ao. O02 O00), 
0.0 0.0 0.0 0.0 0.0 0.0 0.0  OQ2OP" 12 G Mero F050) 0.0 Oe cee cr, 
0.0 0.0 0.0 0.0 0.0 (0.0 “0.0 O00) O07 0 170 70.0) 0. On Ole ter 
0.0 .0.0)..0.0..0.0.,0.0 . O6:0i00 0.10. OS Or On O Om Ol Oe Oe Or Omen 
0.0 0.0 0.0 0.0 0.0 0.0 "0.0 0.0 0.0 ©.0 0.07 = 1. Ge; 04 0.0 
0.0 0.0 0.0 0.0 0-0 0.0 0.0 0.0 0.. 0550.0. 0.0.0. OR Ges 04 aor 
0.0 0.0 0.0 0.0 0.0 0.0 0.0. 0.05 0.0 70.0500, 0500.0 SOO oe... 
0.0 0.0 =.1479€+04 0.0 0.0 =-f1143d+03" 90.0 0-70 9—. 2529da-01) 0-0 

0.0 =.2697d+03 050 0-0; 
0.0 0.0 0.0 -.14914€+04 0.0 0.05 =. 1149d+03"55 0-0 0. 0Ne-. 25560" 0,2e 
0.059.060.6067" =—.27014d+03" -0.0; 

OOP 0.0 O40) 0.0 )0—. 53020705 20-06 .0 Om --1 6 16a..04 O70) -O..0 

~.1790d+02 0.0 0.0 =-.6053d+03] 

bg=(0.0 0.0; 
O20 0.08 
0.0 aor ©, 
0,10 370)..0); 
0:0", <O.c- 0 
OO 2.0. 0+ 
0.0 7500. 
0.0 O26. 
0] 0G. 0: 
070 0.0; 
0.027070 


‘ 
~1479d+04 0.0; 
O2o -1491d+04; 
.5308d+05 0.0] 
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ea-11-0 0.0 0.0 O20 5 900 =-0 1.00. Om 0. O ero ..0/ 80.0 0.0" 0.0" 50.0 0.0; 
ero 1. O50. 0810.0 OF09 0765 0.0550.60 020" 0.0. 0.0 0.0 0.0 0.0] 

dg=zeros(2) 

$pause 

disp(’ se 

Gasp ( ’ Balanced realization of the X-29 state space representation’) 

disp(’ ") 


{agbl,bgbl,cgbl,g,t]=obalreal(ag,bg,cg) 

ag=agbl; bg=bgbl; cg=cgbl; 

$pause 

w=logspace(-3,3,100); 

Svg=sigma(ag,bg,cg,dg,1,w); svg=20*logl0(svg) ; 

semilogx(w,Svg) 

title(’X-29 SV PLOT OF PLANT G’) 

Xxlabel(’FREQUENCY - rad/sec’) 

ylabel(’SV - dB’) 

grid 

meta hx29gsv 

disp(’ 2) 

disp(’ Calculate the poles and transmission zeros of the balanced open’) 
Gasp (’ loop plant) 

disp(’ o 

poleg=eig(ag), tzerog=tzero(ag,bg,cg,dg) 

disp(’ io} 

disp(’ = 

disp(’ Determine determine the condition number of ag’) 

disp(’ E) 

disp (’ ~ 

condag=cond(ag), rcondag=rcond(ag) 

disp(’ ») 

disp(’ oS 

$pause 

disp(’ ‘) 

disp(’ << Design Specifications >> ’”) 

Gisp(’ ’) 

Gusp(’ 1). Robustness Spec. : -40 dB roll-off, -20 db ®@ 100 Rad/Sec.’) 
disp(’ Associated Weighting: ’) 

disp(‘ y) 

disp(’ -1 1000 ’) 

disp(’ W3(Ss) =  ------ * J (fixd) ’) 

disp(’ 2 2x2") 

Gasp { ’ s a) 

disp(’ 9 

Gdisp(’ ’) 

disp(’ 2). Performance Spec.: minimizing the sensitivity function’) 
disp(’ as much as possible.’) 

disp(’ Associated Weighting: ’) 

disp(’ ‘) 

disp(’ a) 
disp(’ -1 =12701(100s + 1) y 
disp(’ W1(s) = Gam * efee-ee---- * TI 
aasp(’ 2X2’) 

disp(’ Gas. + 1) ) 

disp(’ i} 

disp (’ where "Gam" in this design is iteratively updated from 1’) 
w=logspace(-3,3,100); 

k=1000; mn=[2 2]; tau=0.0; 

nuw3i = [0.0 k]; dnw31i = [1.0 0 0}; 


J) 


SVWw31 
nuwli 
svwli 


tout 


bode (nuw3i,dnw3i,w); Svw3i = 20*1l0g10(Ssvw3i); 


[1,0 0.01); dnwil (6.1 120), 
bode (nuwli,dnwli,w); svwli 


20*1l0g10(Svwli) ; 


aw2=-0.025*eye(2) ; bw2=zeros(2); cw2=zeros(2); dw2=-0.025*eye(2) ; 


disp(’ 
disp(’ 
disp(’ 
tpause 


") 
") 


(strike a key to see the plot of the weightings em) 


semilogx(w,svwli,w,Svw3i) 


grid 


title(’X-29 Design Specifications’ ) 
xlabel(’Frequency - rad/sec’) 
ylabel(’1/W1 & 1/W3 - db’) 

text (.01,0,’1/Wi(s)’) 
text(.5,100,71/W3(s)’) 

meta hx29specs 


t pause 
cle 

disp(’ 
disp(’ 
disp(’ 
Gdisp(’ 
disp(’ 
disp(’ 
adisp(’ 
euso. 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
Gispi(? 
aisoi~ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
asp 
also 
adisp(’ 
disp(’ 
disp(’ 
tpause 
ele 

disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
aisp(” 
disp(’ 
disp(’ 
disp(’ 
disp(/’ 
disp(’ 
disp(’ 
disp(’ 


<< Problem Formulation >>’) 


") 
- 


Form an augmented plant P(s) with these two weighting functions: ’) 


1). Gam*W1 penalizing error signal "e"’) 
‘) 


‘ 


2). W3 penalizing plant output "y"’) 


and find a stabilizing controller F(s) such that the Hinf norm’) 


a) 


of TF Tylul is less than or equal to one, i.e.’) 
2 
|Tylul | < or = 1,’) 
F(s) ink) 
where ‘) 

-1 Ae 

Tylul = Gam*W1*(I + GF) = Gam * W1 * S > 
=a W3.* (I2=)S) ie) 


W3*GF*(I + GF) |’) 


-)") 


(strike a key to continue .. 


<< DESIGN PROCEDURE >>) 


ee ke ek ke we ke we wR we we ke ke ke we eR ke eR we eR eR we eR ke ke ke ke wk a) 


{Step 1}. Do plant augmentation (run AUGMENT.M or */) 

AUGX29.M) x!) 
=) 
x’) 
x’) 
x!) 
x!) 
x!) 
Redo the plant augmentation for a *“f) 


higher "Gam" and rerun HINF.M x!) 
i 


{Step 2}. Balance the augmented plant for better 


numerical condition if necessary 
{Step 3}. Do Hinf synthesis with "Gam" = 1 


{Step 4}. 


+ + + + + + + © 4H HB HM FH 
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disp(’ ‘) 

disp(’ (strikesa key to continue wa.) ) 
tpause 
elc 
disp(’ 
disp(’ 
disp(’ 
aisp{’ Assign the cost coefficients "Gam" with Gam=1 ’) 

disp(’ ") 

disp(’ serving as the baseline design ....’) 

Grsp(’ ') 

Q@isp(% rr rrr rrr nnn rrr e- e) 
disp(’ augment % Plant augmentation of the X-29 dynamics’) 
Gisp(’ —— nnn enne ") 
disp(’ 2 

disp(’ ") 

haugx29 
dasp(’ ’ 
adisp(’ ' 
disp(’ (strike a key to continue ...)’) 
$ pause 

crc 

disp(’ “} 

disp(’ 7”) 

$disp(’ Do state space balancing on the augmented plant if needed’) 
tdisp(’ ia) 

tdisp(’ : 

ml abal,bbal ,cbal,g,tj=obalreal(A, [Bl B2},{C1;C2]})} 

Sh=aeal, Bl=bbal(:,1:2), B2=bbal(:,3:4), Cl=cbal(1:4,:), C2=cbal(5:6,:) 


~~ ™~ ~~ 
wee eee ee” 


disp(’ -) 

disp(’ » 

disp(’ The transmission zeros, poles and condition number of the augmented’) 
disp(’ plant follow. In addition, determine if (A,Bl) & (A,B2) are ") 

disp (’ stabilizable and if (C1,A) & (C2,A) are detectable.’) 

disp(’ ) 

disp(’ : 


tzeroaug=tzero(A,[Bl B2)],[C1;C2],{D11 D12;D21 D22)), poleaugA=e1ig(A) 
condaugA=cond(A), rcondaugA=rcond (A) 

eps=eps 

toldef=10*max(size(A))*norm(A,1)*eps 

tol=100*eps*norm([A Bl1}) 

mcr ele, Cle, t,kK)]—ctrbt (A,B1,C1,tol) 

tol=100*eps*norm([A B2)) 

Paee pee, €2c,t,k]=ctrbf (A,B2,C2,tol) 

tol=100*eps*norm([A;C1)) 

[Ate,b1e,Clo,t,k)=obsvf (A,B1,C1,tol) 

tol=100*eps*norm([A;C2])) 

{A20,B20,C20,t,k]=obsvf (A,B2,C2,tol) 

clear condag rcondag poleg tzerog svw3i svwli 

clear condaugA rcondaugA 

Clear functions 

tpause 
disp(’ 
disp(’ 
ee ee = = = = eS = = = -- - _-- -— - = = = = = - a) 
disp(’ Hine $ Running script file HINF.M for Hinf optimization’) 
Ci nn a a we ) 
aretype=’Schur’ 

hinf 

disp(’ ‘) 


i 


disp(’ ”) 


disp(' (strike a key to coneinue 0.) 9 
clear functions 

$pause 

disp(’ ) 

disp(’ -) 

disp(’ State space representation of the full order controller’) 
disp(’ (acp, bcp, ccp, dcp) with its poles and condition number’ ) 
disp(’ o) 


acp, bcp, cep, acp 
polec=eig(acp) 
condacp=cond(acp) 


disp(’ : 
disp(' Minimal realization of the controller’) 
disp(’ e) 


toldef=10*max(size(acp) )*norm(acp,1)*eps 

tol=100*eps*norm({acp bcp;ccp dcp) ) 

(acpm,bcpm,ccpm,dcpm)=minreal(acp,bcp,ccp,dcp) 
") 


disp {' 

disp(’ Balanced realization & model reduction of the minimal controller’) 
disp(’ Ss 

{acpbl,bcpbl,ccpbl,g,tj)=obalreal (acpm, bcpm, ccpm) 

elim=(10} 

(acpr,bcpr,ccpr,dcprj=modred(acpbl,bcpbl,ccpbl1,dcpm,elim) 

disp(' i) 

arvep(/ Poles, controllability, observability, and condition of the ’) 
disp(’ balanced, reduced order controller’) 

disp(’ ny 


poleacpr=eig(acpr) 

tol=100*eps*norm({acpr bcpr}) 

(acpe, bepe, ccpc,t, K)=ctrbf(acpr, bepr,cepr, co 
tol=100*eps*norm((acpr;ccpr])) 

{[acpo, bepo, ccpo,t, Kk ]=onsveéi (acpr, bepr,cepr, tol 


condacpr=cond(acpr), rcondacpr=rcond(acpr) 

acp=acpr; bcp=bcpr; ccp=ccpr; dcp=dcpr; 

disp(’ 2) 

disp(’ CLTF Tylul (acl, bel, cel, dcl) and its poles (reduced order.) 


(acl,bel,ccl,del]=lftf (sysp, dimp, acp,; bep, ceenacr) 
polet=elig(acl) 

$pause 

disp(’ 2) 

disp(’ a) 

hpltopmod % Preparing singular values for plotting 
end 
disp(’ 
disp{’ 
disp(’ (strike a key to continue ...)’) 
$pause 

disp(’ =) 

disp(’ a5) 

disp(' Open loop state space representation of controller/plant series’) 
disp(’ ) 

(algf,blgf,clgf,dlgf)=series(acp,bcp,ccp,dcp,ag,bg,cg,dg) 

polol=eig(algf) 

disp(’ ) 

disp(’ a) 

disp(’ (strike a key to continue ...)’) 
$pause 
disp(’ 
disp(’ 


~ 


eee el” 


= ~ 
eee ee 
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avary h3x29.dat 


format short e 
aisp( ’ n3X29 om”) 
disp(’ ) 
Gisp(’ mpasesecripe file 1s designed toysolve the Hinf optimal control ’) 
aasp(/ problem for the X-29. The 14th order FDLTI model, in state ’) 
disp(’ space form, is that of the X-29 aircraft and actuator dynamics.’) 
disp(’ Two states are those of the aircraft dynamics, i.e., alpha ’) 
aiepi(’ and q. The remaining 12 states are the dynamics of the ’) 
disp (’ three, fourth order actuators, i.e., the canard, flaperon, and ’) 
disp(’ strake actuators. The order of the unbalanced states is as ’) 
@isp(’ follows: ’) 
disp(’ ) 
disp(’ wieha. Gq, de, GL, ds, dcdot, dfdet, dsdot, dedbldot,”) 
Gisp(’ Gtepidaot, adsdpldot, dcetrpldot, dftrpldot, dstrpidot’) 
disp(’ o, 
disp(’ Given the open loop transfer function G(s)=Cinv(IS-A)B+tD, a ’) 
aasp(’ stabilizing controller F(s) will be found such that the Hinf norm’) 
disp(’ of Tylul is less than or equal to one. ‘'") 
disp(’ a) 
disp(’ ” 
pause 
cle 
disp(’ X-29 aircraft and actuator state space representation’) 
disp(’ ‘) 
Gieot’ ") 
ag=[~.4181d+00 .9960d+00 <-.2269d-01 -.1213d+00 -.1948d-01 -.9493d-03 
-4427d-04 Ree oe ~-1451d-05 -.2162d-04 mee Oc > 0-0 O20 0.07 
~-5474da+01 -3424d+00 ,2 90507 O) . =, 1S coc rel -1058d+01 ~-3898d-02 ... 
~.1164d-01 - 6397a- =O 2a O00 = 0 3a = 3620-03 eae Ogu, 08.0 0.0; 
Pn oor Ome. O Omori O Owes O.O  O.0° 9020 0.0 0.0 0.0 90.0; 
oto nO momma 20.0 0.0 2:0 "0.0 0.0 O20, “0.0 0.0 0.0 0.0; 
Poo. O.Om 0.0, 0.0 0.0 0.0 1.0 026 0.0 WakO0 0.0 0.0 0.0; 
fPeemmorO . O-.0° 550.6 0.0 60.0 0.0 0.0 "170" 0.0 0.0 0.0 0.0 0.0; 
Om. Ome O.O0n O2On O20 00> 0,0 5.0.0 1.0 0.0 0.0 0.0 0.0; 
amore” OFC O.0 O.0 0.0 0.0 0.0 650 0.0 2.0 0.0 0.0 0.0; 
remo O.0 10.6 0.0 0.0 0.0 0.60 0.0.0.0 060.0  12:0e+04 0.0 G20; 
faoeeo.o, 0.0 90.0 O.1.0 0.0 0.0 0.0 0.098670 0.0 0.0 1.0e+04 0.90; 
momo 0) 0.0 0.0. 0.0 0.05 0.0 0.0.0.0 0.0 0.0 ae: 0 0.0 1.0e+04; 
0.0 0.0 =-.1479d+04 0.0 0.0 -.1143d+03 0.0 0.0 ~25029G+01 0.0 2.2. 
0.0 .2697d+03 0.0 0.0; 
rowmO.O) 0.0 =.1491d+04 0.0 0.0 =.1149d403 0.0 0.0 -.2536d+01 ... 
peo 0.0 =-.2701d+03 0.0; 
@miowno.O . O.O ~O020 =.5302a€405 0.0 #@©:.0 =-.1816d+04 0.0 0.0 
-.1790d+02 0.0 0.0 -.6053d+03] 
b@=(0.0 0.0 0.0; 
Pao 0.0 0.0; 
O20 0.0 . 0.0; 
OnOm 0. 0 =O, 0; 
mo 0.0 0.0; 
Seon 0.0 0.0; 
Of) 0.0 0.0; 
oon 0.0 0.0; 
ono O.0, O.0; 
oO 80.0 0.0; 
O00 0.0 0.0; 
.1479d+04 0.0 0.0; 
O70) .149i1d+04 0.0; 
0.0 0.0 .5308d+05]; 
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eg={1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 | 02007 0a. Ce oor 
0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 @.0 (GO OU OreG Ooo coo. 
0.0 0.0 1.0 0.0 90.0 0.0 0.0 0.0 0.0 3020 0. 00-000) Cor 

dg=zeros(3,3); 

tpause 

disp(’ ) 

disp(’ Balanced realization of the X-29 state space representation’ ) 

disp(’ 4) 

(agbl,bgbl,cgbl,g,tJ=obalreal(ag,bg,cg) 

ag=agbl, bg=bgbl, cg=cgbl, 

w=logspace(-3,3,100); 

SvVg=sigma(ag,bg,cg,dg,1,w); svg=20*log10(Svg) ; 

semilogx(w,svg) 

title(’X-29 SV PLOT OF PLANT G’) 

Xlabel(’FREQUENCY - rad/sec’) 

ylabel(’SV - dB’) 

grid 

meta h3x29gsv 

disp(’ -) 

disp(’ Calculate the poles and transmission zeros of the balanced open’) 

disp(’ loop plant’) 

disp(’ = 

poleg=eig(ag), tzerog=tzero(ag,bg,cg,dg) 

disp(’ “ 

disp(’ a) 

disp Determine determine the condition number of ag’) 

disp (’ ) 

disp(’ ‘:) 

condag=cond(ag), rcondag=rcond (ag) 

disp(’ 1) 

disp(’ ‘) 

pause 

Grsp (os) 

disp(’ << Design Specifications >> "’) 

disp(’ ’) 

Gisp(" 1). Robustness Spec. : -40 GB roll-off, -20 db @ 100 Rad/Sec.’) 

disp(’ Associated Weighting: ’) 

disp (* “) 

disp(’ -1 1000 ’) 

disp(’ W3(S) = cecc-- * T eiixd)) 

disp(’ 2 3x3’) 

disp(’ Ss 3) 

disp (’ ‘) 

disp(’ ’) 

disp(’ 2). Performance Spec.: minimizing the sensitivity function’) 

disp(’ as much as possible.’) 

disp(’ Associated Weighting: ’) 

aasp(* 7) 

disp(’ oS) 

disp(’ =) =“) .01(100s" +" 1) *) 

disp(’ W1(s) = Gam * errrcres---= * I’) 

Gisp(’ 3x3‘) 

disp(’ (.1s + 1) “) 

disp(’ ee, 

disp(’ where "Gam" in this design is iteratively updated from 1’) 


w=logspace(-3,3,100); 
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k=1000 
nuw3i 
SVw3i 
nuwli 
Svwli 
aw2=-0 
disp(’ 
disp (’ 
disp(’ 
tpause 
semilo 
grid 
title( 
xlabel 
ylabel 
text (. 
text(. 
tmeta 
tpause 
elc 
disp(’ 
ayen(’ 
asp (’ 
dasp( ’ 
disp(’ 
disp(’ 
disp(’ 
dise:( ’ 
aisp( ’ 
disp (’ 
disp(’ 
aso ’ 
disp (’ 
disp(’ 
disp(’ 
disp’ 
aarsei( ’ 
disp(’ 
pause 
Guc 
asp (’ 
airs (’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 
adres (’ 
disp(’ 
disp(’ 
eis” 
disp’ 
disp(’ 
disp’ 
disp(’ 
disp(’ 
anspit ’ 
ears ’ 
pause 
cic 


mn=(3 3}; tau=0.0; 

[0 O kK}; dnw3i [2530 SOON; 
bode (nuw31,dnw3i,w); svw3i 
ple 0.0.01); dnwil =(0.171.01; 

bode (nuwli,dnwli,w); svwli = 20*1l0g10(svwli) ; 
~-O25*eye(3); bw2=zeros(3); cw2=zeros(3); dw2=-0.025*eye(3) ; 
c 

| 


a 
= 


20*10g10(Svw31) ; 


iow Ww Ws. 


(strike a key to see the plot of the weightings sa) 


gx(W,SVWli,w,SvVw3i) 


’X-29 Design Specifications’) 
(‘Frequency - rad/sec’) 
(/1/W1 & 1/W3 - db’) 
O10, Wy Wis)” ) 

>, 100, 1/W3(s)") 

h3x29spc 


; << Problem Formulation >>’) 

Form an augmented plant P(s) with these two weighting functions: ’) 
1). Gam*W1 penalizing error signal “e"’) 

2) nie mpenallzing control input “yl” ) 

2). W3 penalizing plant output "y"’) 


and find a stabilizing controller F(s) such that the Hinf norm’) 


of TF Tylul is less than or equal to one, i.e.’) 
‘) 
|Ty1u1| <or-= 1," ) 
F(s) rit) 


-)’) 


(strike a key to continue 


<< DESIGN PROCEDURE >>’) 


ze ek ke ek eR ke ke eR eR ke Re ke ke RR ke eR eR eR ke ke Rk eR ke ke KRY 
* [Step 1}. Do plant augmentation (run H3AUGX29.M ) x!) 
* *!) 
* [Step 2}. Balance the augmented plant for better */) 
* numerical condition if necessary *?) 
* */) 
* [Step 3}. Do Hinf synthesis with "Gam" = 1 =e) 
* a!) 
* [Step 4]. Redo the plant augmentation for a — 
* higher “Gam" and rerun HINF.M a 
kk ke ke ke ke ke ke ke ke eR Rk ke Rk eR Rk eR Rk Rk eR Rk eR ke RY) 


=~ ™ 


a eee” 


-)") 


(strike a key to continue .. 
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disp( 
disp(’ 
disp( 
Gisp(’ Assign the cost coefficients "Gam" with Gam=1 ') 

disp(' i 

disp(’ serving as the baseline design ....’) 

arsp( a) 

Gisp($ —_ wn rem re rn nn ee ee ee a) 
aispie augment %* Plant augmentation of the X-29 dynamics’) 
Gisp($ «—_—-_—_ mtn mn eer nm err nn en nn nnn nnn ee ee - os) 
disp(’ 
disp(’ 
h3augx29 
disp) 

disp(’ ‘) 

disp(’ (striké a key to continue =...) ) 
tpause 

clic 

disp(’ a) 

disp(’ ) 

$disp(' Do state space balancing on the augmented plant if needed’) 
tdisp(’ ) 

$disp(’ v4 

%$(abal,bbal,cbal,g,tj=obalreal(A,[Bl B2],({C1;C2]}) 

tA=abal, Bli=bbal(:,1:2), B2=bbal(:,3:4), Cl=cbal(1:4,:), C2=cbal(5:6,:) 


= SS & 


- 


ee 


disp(' “) 

disp(’ <) 

disp(’ The transmission zeros, poles and condition number of the augmented’) 
disp(’ plant follow. In addition, determine if (A,B1) & (A,B2) are ’) 
disp(’ stabilizable and if (C1,A) & (C2,A) are detectable.’) 

disp(’ -) 

disp(’ : 


tzeroaug=tzero(A, [Bl B2],[C1;C2],[(D11 D12;D21 D22}), poleaugA=eig(A) 
condaugA=cond(A), rcondaugA=rcond (A) 

eps=eps 

toldef=10*max(size(A))*norm(A,1)*eps 

tol=100*eps*norm([A Bl1])) 

[Aic,Blc,Cic, t7x)—ctrcbf (A, Bl ei cer) 

tol=100*eps*norm((A B2}) 

[Ade, B2c;C2c, tiaeetrbt (A, Ba, eentoly 

tol=100*eps*norm([A;C1}) 

[Aio,Blo,C1lo,t,k)]=cebsvi(A, Ei eice.) 

tol=100*eps*norm([A;C2]) 

[A2Zeo, B20,C20,t,K]=obsvti (A, B22) c2, co, 

clear condag rcondag poleg tzerog svw3i svwli 

clear condaugA rcondaugA 

clear functions 

tpause 
disp(’ : 
disp(’ 7 
Qisp($ rt rr SS Se ee Oe ee a ee oe ee ) 
adisp(’ hae ¢ Running script file HINF.M for Hinf optimization’ ) 
Gisp($ rr 8 8 8 = = = = = =) 
aretype=’Schur’ 

hing 

atep(’ ") 

disp(’ ‘) 

disp(’ (strike a key to continue ...)’) 
clear functions 

$pause 
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asp (’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 


a2 
©, 

State space representation of the full order controller’) 
(acp, bep, ccp, dcp) with its poles and condition number’ ) 
’) 


a@ep,bcp,ccp,dcp 
polec=eig(acp) 
condacp=cond (acp) 


disp { ’ 
a@asp(’ 
disp(’ 


2) 
Minimal realization of the controller’) 
a 


toldef=10*max(size(acp) ) *norm(acp,1)*eps 
tol=100*eps*norm([{acp bep;ccp dcp}) 
(acpm, bepm, ccpm,dcpm)=minreal(acp,bcp,ccp, dcp) 


disp(’ 
aqisp(’ 
disp(’ 


c 


Balanced realization & model reduction of the minimal controller’ ) 


é 


(acpbl,bcpbl,ccpbl,g,t]j=obalreal(acpm, bcpm, ccpm) 


elim=([10}] 


Paepu, bepr, ccpr,dcpr ])=modred(acpbl,bepbl,ccpbl,dcpm,elim) 


aiso(’ 
aisp(’ 
disp(’ 
Gisp(’ 


c 


Poles, controllability, observability, and condition of the ') 
balanced, reduced order controller’) 


c 


poleacpr=eig(acpr) 

tol=100*eps*norm([{acpr bcpr]) 

meme, Ocpe,ccpc,t,kK}]=ctrbf(acpr,bepr,ccpr,tol) 
tol=100*eps*norm({acpr;ccpr})) 

feero, -epo,ccpo,t,ki=obsvf(acpr,bcpr,ccpr,tol) 


condacpr=cond(acpr), 
acp=acpr; 


disp(’ 
disp(’ 


rcondacpr=rcond (acpr) 

bep—bepr; cep=ccpr; acp=<dcpr; 

o) 

Sith elviuieracl bel, eccl,-del). and ats poles (reduced order)’) 


feeiebecl, col, dclj|=lftf(sysp,dimp,acp,bep, cep, dcp) 
polet=eig(acl) 


pause 
disp(’ 
disp(’ 
neplot 
end 

disp(’ 
disp(’ 
disp(’ 
$pause 
disp(’ 
arep(’ 
disp(’ 
disp(’ 


‘) 
*) 


% Preparing singular values for plotting 


~» ~ 


(strike a key to continue ...)’) 


i) 
‘) 


Open loop state space representation of controller/plant series’) 
é 


feat, blgt,clgf,dlgf j=series(acp, bcp,ccp,dcp,ag,bg,cg,dg) 
polol=eig(algf) 


disp(’ 
disp(’ 
G1sp(’ 
pause 
aise’ 
disp(’ 
aisp(’ 
ais (’ 
disp(’ 
disp(’ 


") 
+) 
(strike a key to continue ...)’) 


) 


Closed loop state space representation of controller/plant series,’) 
controllability, observability, and condition number of the closed’) 


loop acgf matrix. ’) 


) 


63 


disp(’ =) 

[acgf ,begf , cegt ,dcegFf |=feedbk (alagt bilagt cigt ,digr,. 
tol=100*eps*norm([acgf bcgf)) 
[acgfc,bcgqtic, ccafce, ©, k)—-ctrbr (acer, pcgtf ,ccgt tom 
tol=100*eps*norm( [acgf;ccgf)) 

facgfo, begfo, ccaro, ©, kK | =enbsvVi(aeqtbegt,ccat ea) 
condacgf=cond(acgf) 


disp(’ i) 

disp(’ &) 

disp(’ (strike a key to continue ...)’) 
pause 

disp(’ ) 

disp(’ Poles of the closed loop system’) 
disp(’ ‘) 

polcl=eig(acgf) 

h3aly 

th3rsp 

end 


Plant Augmentation for the X-29 Hinf problem as W3 is not a 

Peeper transter function. Includes contingency for adding W2 to 
Ensure di2 is full column rank. This script file is designed for 
the X-29 system without theta as a state, ie, 3 inputs & 3 outputs. 
tdisp(' ry 


of OF OF OF OF 


tdisp(’ z) 

adisp(’ << Plant Augmentation >>’) 
$Gam=gama(1,i) 

Gam = input(’ Input the cost coefficient "Gam" = ”) 
cgb = 1/k*[(cg(1,:) *ag*ag; 


cg(2,:)*agtagtag*taut+cg(2,:) *ag*ag; 
cg(3,:)*agtagtag*tautcg(3,:) *ag*ag); 
Pee egq(l,:)*ag*bg; 
cg(2,:) *agt*tag*bg*tau; 
cg(3,:) *agtag*bg*tau); 
Reve—eGan*|anwl1l;0 0;0 0;0 O-dnw1i;0 0;0 0:0 O;dnwli}j; 
dwl = nuwli; 
SySw2=[aw2 bw2;cw2 dw2); xw2=3; 
Meee be, cl,cz2,D11,D12,D21,D22 )=haugmod(ag,bg,¢g,cgb,dg,dgb,nwl1,dwl, mn) 
fae biz, cl,¢2,D11,012,D21,D22 }=haugx29pl1 (ag, bg,cg,cgb,dg,dgb,nwl1 ,dwl,sysw2, xw2, 


dgb 


ecdispie *) 

tdisp(’ - - - State-Space (A,B1,B2,C1,C2,D11,D12,D21,D22) is ready for’) 
tdisp(’ the Small-Gain problem - - -’) 

% 

$ ------ End of AUGX29.M ----- %$ *Z 
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disp(’ 
disp(’ 
disp(’ 
disp(’ 
disp(’ 


Closed loop state space representation of controller/plant series,’) 
controllability, observability, and condition number of the closed’) 
loop ecGi matrix...) 


uf 


facgf, begf,ccgf,decf |—feedbk (alot tigt claf dict 2) 
tol=100*eps*norm({acgf bcgf}) 
facgfc,bcgfic,;cegte,t,k |=etrbt (acgi,ceat cegrro 
tol=100*eps*norm([(acgf;ccgf]}) 

{aegto, bcgfo, cegfo,t,k ]—obsvi (aegt, begt, cegr 1.6) 
condacgf=cond (acgf) 


disp(’ 
disp(’ 
disp(’ 
t pause 


‘) 
") 


(strike a key to continue ...)’) 


tdiary hx29.dat 


disp (” 
disp (~ 
disp(’ 


c 


Poles of the closed loop system’) 
c 


polcl=eig(acgf) 


thx29aly 
thx29trsp 


end 
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format short e 

w=logspace(-3,3,100); 
(magl1,phasel]=bode(algf,blgf,clgf,dlgf,1,w); magl=20*1l0g10(magl) ; 
{mag2,phase2]=bode(algf,blgf,clgf,dlgf,2,w); mag2=20*1l0g10(mag2) ; 
[mag3,phase3]=bode(algf,blgf,clgf,dlgf,3,w); mag3=20*1l0g10(mag3) ; 
semilogx(w,magl(:,1),w,mag2(:,1),w,mag3(:,1)) 

title(’X-29 OPEN LOOP (GF) BODE PLOT OF " alpha " ”) 
xlabel(’FREQUENCY - rad/sec’) 

ylabel(’GAIN - dB’) 


grid 

$meta h3bdl 

$pause 

semilogx(w,magl ( eve, MAGe. «5 2), 6, 0ags(: ,2))) 
title (’ X—-29 OPEN LOOP (GF) BODE PLOT OF “q " *) 


xlabel(’FREQUENCY - rad/sec’) 

ylabel(’GAIN - dB’) 

grid 

tmeta h3bd2 

$pause 

[(mag1, phaselj=bode(acgf,bcgf,ccgf,dcgf,1,w); magl=20*1log10(magl)j; 
[mag2, phase2 ]}=bode(acgf,bcgf,ccgf,dcgf,2,w); mag2=20*log10(mag2) ; 
(mag3, phase3 )=bode(acgf,bcgf,ccgf,dcgf,3,w); mag3=20*1log10(mag3) ; 
semilogx(w,magl(:,1),w,mag2(:,1),w,mag3(:,1)) 

title(’xX-29 CLOSED LOOP BODE PLOT OF " alpha " ”) 
xlabel(’FREQUENCY - rad/sec’) 

ylabel(’GAIN - dB’) 


grid 

text(15,-20, ’-r--- TO. ri!) 
text (.22,-20,'----- TO £2’) 
text(.1,-60, ’----= TOES’ ) 
meta h3bd3 

pause 


semilogx(w,magl(:,2),w,mag2(:,2),w,mag3(:,2)) 
title(’X-29 CLOSED LOOP BODE PLOT OF "gq "* ? 
xlabel(’FREQUENCY - rad/sec’) 

ylabel(’GAIN - dB’) 


) 


grid 

text (.1,-10,'’----- LOmr2") 
text(.1,730, '-r--- TO@E 1") 
text(.1,°58,'’----- LOTT”) 
meta h3bd4 

$pause 

disp(’ ‘) 

end 
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APPENDIX B 


State Space Realization of the Open Loop Plant G 
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APPENDIX C 


State Space Realization of Reduced-Order Controller for 2 /O System 


acpr = 


Columns 1 through 6 


-1.0399e-02 
-1.1271e-04 
-1.0340e-01 
-5.7615e-02 

7.8601e-02 
-3.4048e-03 


-4.9864e-03 
3..65939e-03 
4.7908e-02 
9.2037e-03 
8.7200e-04 
8.3046e-04 
8.2445e-04 
1.7763e-04 
2.3188e-05 


-1.2262e-05 
=. olL96eE—02 
-2.9977e-02 
6.0104e-02 
1.3419e-02 
-3.6138e-02 


-3.9735e-02 
-2.9112e-02 
-1.1300e-02 
-1.4652e-03 
6.8844e-03 
-3.7921e-03 
1.0942e-04 
9.0089e-04 
3.6527e-05 


Columns 7 through 12 


2.5870e-03 
-4.0336e-02 
9.2118e-01 
1.0142e+01 
-5.1779e+01 
3.2298e+00 
-1.3998e+01 
=3.70)2e+01 
2.0492e+01 
1.5648e+00 
4.4792e+00 
-2.3356e+00 
2.7969e-01 
5§.6035e-01 
3.03 10e=02 


9.1286e-03 
2.5238e-02 
3.4965e+00 
-1.3763e+01 
4.6467e+01 
1.4770e+01 
3.5254e+01 
-1.4396e+01 
-7.6498e+01 
-4.6369e+00 
6.6710e+00 
-3.8924e+00 
-2.1520e-01 
9.1670e-01 
2.5482e-02 


Columns 13 through 15 


-4.8787e-04 
-1.1798e-03 
-1.8641le-01 
6.2164e-01 
-1.9086e+00 
-1.9215e-01 
-7.3742e-01 
1.4035e+00 
-~5.7313e+00 
2.3404e+00 
-8.1543e+00 
2.4647e+00 
-4.5631e+01 
3.3396e+01 
~2.7014e+00 


-~4.7600e-04 
“4.6223e-04 
-1.7785e-01 
3.9705e-01 
-1.5949e+00 
9.8820e-02 
-2.4886e-01 
8.0952e-01 
-~1.7434e+00 
3.5305e+00 
-2.5954e+00 
1.5106e+01 
-6.3440e+01 
-7.1315e+01 
-3.0793e+00 


-1.0460e-01 
-4.6795e-03 
-1.9525e+01 
-2.7029e+01 

3.1746e+01 
-2.7492e+00 


-3.1302e+00 
2.9037e-01 
1.701 0e+01 
3.2877e+00 

-9.9621le-02 
1.8196e-01 
3.0302e-01 

-3.6094e-02 
9.5788e-03 


8.3886e-03 
~8.4033e-02 
2.3161e+00 
1.9246e+01 
1.9853e+01 
-4.2332e+01 
-6.3778e+01 
9.7679e+01 
-2.4305e+02 
-8.7073e+01 
1.6807e+01 
-1.1275e+01 
-8.1078e+00 
2.5024e+00 
-2.0105e-01 


-1.1648e-05 
-5.0636e-05 
-4.5815e-03 

2.1676e-02 
-5.4329e-02 
-1.3922e-02 
-3.2950e-02 

5.1750e-02 
-2.6308e-01 

1.1476e-02 
-3.65l1le-01 
-1.0758e+00 
-3.3130e+00 
-6.2957e+00 
-1.4418e+02 
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FPNAWHrHW ONWNN OD 


1 
AWN 


=) 5 


=3. 


=e 
= 
=3 
-2. 
4. 


=O 
-8. 


-2. 


-8094e-02 
-3252e-02 
-6394e+01 
-5674e+01 
-6483e+01 
-9997e+00 


-5764e+00 
-0513e+01 
.-2813e+01 


3758e+00 


-1329e+00 
-4002e+00 


1368e-01 


#198 1e=—071 
-8832e-03 


-6940e-03 


5085e-03 


-1564e+00 


0766e+00 


-4744e+01 


6418e+00 
0628e+00 
0219e+00 
9439e+01 
5600e+01 


-0142e+01 


7936e+00 
0534e+00 


-2396e+00 


S73 1e=07) 


=/.6875e-02 
=-3.1861e-02 
=-3.1280e+01 

4.0396e+01 
-1.2030e+02 
~2e9o7z2erTOl 


5§.9253e+01 
-2.4717e+01 
-1.4766e+02 
-2.8661le+01 
1.6446e+00 
-1.9663e+00 
-2.5218e+00 
4.0668e-01 
-7.5896e-02 


5§.3095e-04 
-7.0034e-03 
1.5463e-01 
1.8510e+00 
-1.1656e+00 
-3.2832e+00 
-4.9278e+00 
5.0051e+00 
-4.1233e+01 
-2.4619e+01 
-2.7237e+01 
3.2961e+01 
1.9334e+00 
-7.3070e+00 
-2.0302e-01 


1.3274e-02 
-2.8000e-02 
4.9271e+00 
-7.7033e-02 
5.3509e+01 
-9.8984e+00 


-2.4777e+01 
-2.1685e+01 
1.0290e+01 
8.9049e-01 
3.5846e+00 
-1.8969e+00 
1.7864e-01 
4.5402e-01 
2.2842e-02 


in OOSe—0S 
-2.6531e-03 
6.2854e-01 
-1.3642e-01 
4.1721e+00 
-2.1705e+00 
-2.0558e+00 
5.6593e-01 
-1.7944e+01 
~2.1083e+01 
-2.5895e+01 
-1.2550e+02 
-2.7568e+01 
5.7380e+01 
1.1700e+00 


bepr = 


5.2502e+00 
-3.2830e-01 
2.5970e+01 
1.5589e+01 
-1.9609e+01 
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Closed Loop State Space Realization of Controller/Plant Series 
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